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Abstract. In this paper three different approaches towards a 
continuum theory of fiber suspensions are discussed. The first one 
is the classical Thermodynamics of Irreversible Processes with in- 
ternal variables. It derives constitutive equations for fiber suspen- 
sions on the basis of ONSAGERs phenomenological coefficients, 
which are related to the mechanical properties of the fibers. Sec- 
ondly another method of exploiting the dissipation inequality, the 
method introduced by LIU is applied in order to derive results on 
constitutive equations. Finally a mesoscopic background theory is 
discussed. In this approach a distribution of different fiber orienta- 
tions and fiber deformations is assumed. The fiber orientation and 
deformation are additional variables in the domain of field quanti- 
ties. The new field quantities on the enlarged set of variables obey 
balance equations. The mesoscopic balance of mass results in an 
equation of motion for the distribution of fiber orientations and 
deformations. The usual macroscopic fields of continuum mechan- 
ics are averages of mesoscopic fields over the different orientations 
and deformations. Macroscopic quantities characterizing the ori- 
entational order, the alignment tensors are defined by the help of 
the orientation distribution function. They are internal variables 
in the macroscopic theory characterizing the distribution of fiber 
orientations. 

The result of the Irreversible Thermodynamics approach mainly 
concerns the stress tensor. Cauchy's stress is usually assumed to be 
symmetric, but in fiber suspensions it has an anti-symmetric part, 
too. The skew-symmetric part of the stress shows up in the balance 
equation of angular momentum; it is related to the torque density, 
the rate of internal moment of momentum, and the couple stress. 
The latter is a result of the internal structure of the material. The 
fibers are assumed to be microscopic and the stress field smooth 
on macroscopic length scale. If the local structure of the flow- 
ing medium may be characterized by second order pseudo (axial) 
tensors as internal variables then in the linear Onsager equations 
there appears couple stress coupled to the gradient of the angular 
velocity. 

The results of the Liu-procedure are consistent with this obser- 
vation. The couple stress can be calculated as a combination of 
derivatives of the free energy density (see equation 1)3.31 [0 . which 
may have an anti-symmetric part, too. This result may be used 
to calculate the couple stress from an ansatz for the free energy 
density, taking into account the deformation of fibers. 

Another result of the Liu-exploitation of the dissipation inequal- 
ity is that the assumption of the entropy flux being heat flux over 
temperature does not contradict the second law of thermodynam- 
ics, but it does not necessarily follow from the dissipation inequal- 
ity. 



1. Introduction 



A phenomenological theory of suspensions and emulsions of weakly 
deformable spheres has been presented in pQ. The equations of mo- 
tion for the particle deformation, and the constitutive equation for the 
stress tensor have been derived in the whole range of concentrations of 
the emulsion. They are generalizations of earlier work on dilute suspen- 
sions E] and emulsions [U |3] . Thermodynamic arguments have been 
applied in order to reduce the number of unknown constitutive coeffi- 
cients. In semi-dilute suspensions hydrodynamic interactions between 
particles are important and their contribution to the stress tensor has 
been investigated. As it is expected for all materials with an internal 
structure the resulting stress tensor has an anti-symmetric part. An 
overview over calculations of the stress tensor of suspensions of elas- 
tic fibers of fixed orientation as well as on the stress tensor of rigid 
orientable fibers can be found in j^j. 

A constitutive model taking into account the possible damage of 
fibers of a glass- fiber woven polyester composite is presented in [7j. 
The material is viscoelastic with a dependency of elastic properties on 
the micro-structural damage. The constitutive model is derived within 
the framework of continuum mechanics and thermodynamics with an 
internal variable. The phenomenological constitutive coefficients have 
been obtained from experimental tests. Other experimental applica- 
tions can be found in [HI El UH] • 

The derivation of constitutive equations on the ground of kinetic 
theory is similar for polymers (see f. i. ^T]) and suspensions of fibers. 
An application of the stress tensor derived by Batchelor [T21 IT3] from 
a kinetic background theory to corner flow of suspensions of orientable 
rigid rods can be found in [T3] . 

Finally the methods applied in the present paper, namely contin- 
uum mechanics with internal variables and Thermodynamics of Irre- 
versible Processes has been applied successfully in other fields: liquid 
crystals [T3], plasticity of metals and other solids [HI], polymer solu- 
tions. Composites with inextensible fibers have been treated in 
with continuum thermodynamics, exploiting the dissipation inequality 
according to Liu jTSJ. In contrast to the present paper, where the same 
method is applied in the second part, in ^7] the material is a solid fiber 
composite. 

In the last section the methods of the so called mesoscopic theory 
are applied. This theory has been developed for liquid crystals [THl I2D1 
1221- For an overview over other possible applications see 124j . 
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2. Irreversible Thermodynamics of flexible fibers 



The thermodynamic description of systems with dynamic (internal) 
variables has been applied with success to several kind of physical phe- 
nomena such as electric conduction [25] electric and magnetic polariza- 
tion [2U 123 I2E] heat conduction and radiation [21 EDI El E2 HUH US] 
viscoelastic and plastic deformations |33 EZl EHl E3 EOJ HU H2J HSJ 
EH US! EH! Ell EEl HI, liquid crystals jSO]. The thermodynamic the- 
ory of visco-elasticity is based on dynamic variables which are second 
order symmetric tensors. They are related to the deformation of elas- 
tic micro-particles or to the distribution function of the orientation of 
rod like or disk like rigid particles. For example some proteins, certain 
saccharoides, and liquid crystals are well approximated by the above 
model, but it fails, if particles are like elastic slender bars, e.g., in the 
case of fiber suspensions, in which the particles may be twisted and 
bent; their affine deformation is not relevant. 

In the first part the flow induced anisotropy is neglected, and the 
orientation distribution of the fibers is presumed isotropic. The purpose 
here is to investigate the consequences of the flexibility of the particles. 

The deformation of a particle (bend and twist) stores energy, so 
influencing the macroscopic mechanical properties of the fluid. 

In the following we will denote by the symbol : the contraction of 
second order tensors over all indices, in components: A : B = A^B^. 
The scalar product between vectors a and b is denoted by a ■ b. 

2.1. Deformation of a Fiber. In our model the fibers are assumed 
to be straight if not loaded. Then one can chose a coordinate s along 
this fiber orientation and an orthogonal tensor U(s) describing the dis- 
tortion of the fiber. The macroscopic angular distortion tensor defined 
by 

(2.1) V ^U-.f 

as 

takes into account the local deformation of the flexible fibers, and is 
the average over all fiber orientations in the continuum element, s 
is the local coordinate along the fiber, and x is the position of the 
continuum element, s is only introduced in order to describe the local 
fiber deformation. 

The tensor ip is obviously skew-symmetric as U is orthogonal, 

(2.2) U T -U = S, 
where 6 is the unit tensor and 

, , (1U T rp dU 

i.e. 



This way, we introduce the angular distortion vector (the vector 
invariant of the angular distortion tensor) as 

(2.5) <p x 6 = cp. 

Only in this case we will denote the vector by the symbol " in order 
to distinguish it from the tensor cp. 

Let n denote the unit vector tangential to the undeformed fiber. The 
scalar product of <p and n gives the twist: 

(2.6) t = 0-n, 

and the component of perpendicular to n is the bend: 

(2.7) h = 0-u{u-0) 

As the fibers are elastic, two material coefficients fi t and fit have to 
be introduced to describe their stiffness. The torque acting in a cross 
section of the fiber is given by 

(2.8) r = fi t n(n ■ 0) + fi b [0- n(n • 0)\ 

The coefficient \i t is the torsion stiffness, and /if, is the bend stiffness of 
the fiber. The elastic energy stored per unit length of a fiber is given 
by 

(2.9) u e = ^(0.nf + ^-(0.nn 

This way, we have obtained the expression for the energy stored by a 
deformed fiber. 

The fibers are deformed by flow and they continuously relax. The 
angular distortion of an individual fiber depends on the local orien- 
tation of the fiber, 

(2.10) 0=0(n). 

The function is odd, (p{— n) = —0(n), according to the definition of <p 
eq. (j2.1|) . as turning the fiber (changing 0(n) to 0(— n)) gives a - sign 
in the line element ds. 

An obvious approximation for the function 0(n) is the first term of 
its expansion by spherical harmonics; 

(2.11) <p(n)=a-n, 

where a is a second order pseudo tensor as is an axial vector. VER- 
HAS [in] elaborated the theory based on this linear approximation and 
the theory gives only one relaxation time, while the relaxation times 
for twist and bend are expected to differ. As it will be shown in this 
section, the above shortcoming is ceased if the approximation (|2.11j) is 
improved with a third order term (with a traceless tensor (3): 

(2.12) 0(n) = a ■ n + n(n • f3 ■ n), 
which can be decomposed as 

(2.13) 0(n) = a°n + ct a ■ n + a. d ■ n - n(n • ct d ■ n) + n(n ■ T ■ n), 

where a = 1/3 tr a, a a is the skew symmetric part of a, ct d is the 
deviatoric part of ol and T is the deviatoric part of the sum of the 



the first and the last term are due to twist, while the others are due to 
bend. 

Inserting the approximation equation (|2.12|) . the elastic energy of a 
fiber per unit fiber length reads 

(2.14) u e (n) = ^(a° + n-T-n) 2 + tj(a a -n + a d -n-n(n-a d -n)) 2 , 
the average of which over all possible fiber orientations is 

u em = ~r f u e (n)d 2 n = 
4vr J S 2 

(2 , 5) = | (> + JL T , T ) + £ (l a > , Q - + ^ , a * 

if uniform distribution of the fiber orientations is supposed, i.e. the 
fiber orientations are distributed isotropically. A more refined descrip- 
tion introducing a non-uniform distribution of orientations is discussed 
in the third chapter. 

If the total length of the fibers in a unit volume is denoted by Z, we 
obtain 

(2.16) u nd = lu em 

for the part of the specific internal energy stored in the deformation 
of the fibers, and not being dissipated. From the formula, one can see 
that this not dissipated part of the energy is a homogeneous quadratic 
function of the scalar a , and the tensors at a , ot d ) and T. 

2.2. Balance equations. The balance equation for mass is 

(2.17) ^ + p Vv = 0, 

which, for volume preserving motions (incompressible liquids) , reduces 
to 

(2.18) V-v = 0. 

The derivative 4 is the material (or total) time derivative 

(2.19) i = i + ^ 

The balance equation of linear momentum is the usual one 

(2.20) p^ = pf + V-t 

at 

where t is CAUCHYs stress tensor, f is the body force per unit mass, 
p is the density, and v is the velocity of the material. The definition 
for the divergence used here reads in components: 

(2.21) (y-t) i = v j t ij . 

Although the internal angular momentum and the torque will be 
neglected later, we write down the balance for the angular momentum 
in general [221 EH] • 

ds 

(2.22) o— = 2w(t) + om + V • II 



Here s is the internal angular momentum per unit mass, m stands for 
the torque exerted on the material by external fields and II is the cou- 
ple stress tensor. The vector w(t) stands for the vector invariant of 
Cauchy's stress tensor. In the following we suppose that the average 
rotation of the fibers (in addition to thermal fluctuations) is generated 
by the macroscopic mechanical motion. Consequently, the internal an- 
gular momentum (due to thermal motion) can be neglected compared 
to the moment of momentum (due to macroscopic motion). Moreover, 
the rotational energy of internal angular momentum need not be re- 
garded. If we also neglect the torque, which is reasonable as long as no 
electro-magnetic fields are present, equation (|2.22|) reduces to 

(2.23) 2w(t)+V-II = 

The balance equation of internal energy is of basic importance for con- 
stitutive modelling. For the total energy, we can write 

^- I p(u + — )dV = - (f q - dA+ (£v-t-dA + 



dt J rv 2 

(2.24) + j uj n -U-dA + J pvfdV + J pu m -mdV 

The left hand side is the change of the total (internal plus kinetic) 
energy of the part of the material regarded. Kinetic energy due to 
particle rotations has already been neglected. The terms on the right 
hand side are heat flow, and power of stress, of couple stress, of forces, 
and of couple forces. The vectors Wn and u> m are the angular velocities 
of what the couple stress and the torque act on, respectively. 

Now, we assume that both u) m and o; n are equal to the macroscopic 
angular velocity of the body, 

(2.25) u) m = wn = uj = -V x v. 

Making use of the other balance equations and of Gauss' theorem, we 
get 

du ° 

(2.26) p— + Vq=d:t d + II:Vu> 

( I 1 

where d and t stand for the symmetric part of the velocity gradient 
and of Cauchy's stress tensor, respectively. 



2.3. Entropy. The independent state variables in the domain of the 
entropy are the internal energy, the pseudo scalar a and the tensors ct a , 
cx. d and T; the latter taking into account the state of the deformation 
of the fibers. 

Assume that the entropy depends on the dissipated energy only, 



(2.27) 



r] = n e (u u nd ). 



where u n d is the non-dissipated energy (stored in the deformation of 
the fibers). From here and from equations (|2.15j) and (J2.16j) we get 
(2.28) 



7] = 7] e 



(u - — —T : T] - l ^\-cx a : ot a + -cx d : ct d ]\ . 

\ 2p l 15 1 2p l 3 5 V 

The balance equation for the entropy, for vanishing entropy supply, 



is 



(2.29) P^ + V- = ct>O 

The entropy flux is supposed to obey the constitutive equation for 
local equilibrium systems 

1 1 ds 

(2.30) d> = —q where — = -— . 

Making use of the entropy function in equation ()2.28|) and the balance 
equation of the internal energy, we obtain the expression of the entropy 
production density 

a = -\d:t d + U:Vu:-lfi t a a -^T:f- 
T{ P 15 

(2.31) --^ot a : q« - : a d - -q ■ VT 
v ; 3 5 T H 

Here T, cx a , and ct d are the co-rotational time derivatives of the ten- 
sors T, ct a , and ot d , respectively. The introduction of co-rotational 
derivatives is here arbitrary, because it holds 

(2.32) T :T = T : T 

(2.33) a° : d a = a a : ct a 

(2.34) ct d : ot d = a d : ct d . 

It is motivated by the fact that later fluxes and forces should be ex- 
pressed in terms of co-rotational time derivatives, i.e. the time deriva- 
tives in a co-moving frame: 

= dT m m „ da a o da d d 

T = —+T-UJ-UJ-T, ct a = — — hex a -uj-uj-ct a , ct d = — — hex d -uj-uj-a c 

at at at 

The tensors II and Vw are decomposed into their symmetric and 
anti-symmetric parts: 

(2.35) n:Vw = tf:(Vw) d + n n :(Va;) a . 

It results the entropy production density 

Ta s = d : t d + U d : (Vw) d + n a : (Vw) fl - l/i t a°a 

(2.36) -^T : f - l -^cx a : a« - : a* - q • — VT. 

v ; 15 3 5 T 

The first term on the right hand side refers to affine deformations, the 



of the distortions of the fibers while the last term is due to heat prop- 
agation. The set of process rates is 

{q; t d ; U d , a<*; T; IT; a»; a } 

and the set of corresponding forces is 

| - ^VT; d; (V«) ; — —a. ; ~j^ T , 

(Vu;)*; -Z^a°; } 

2.4. Constitutive Equations. In an isotropic material, with the above 
thermodynamic fluxes and forces ONSAGERs linear laws read ( |54[ I55| 

H3 EH EB1 EH HIJ HED E2 ) 

(2.37) q = —i^L q VT = -AVT, 

(2.38) t d = 27?<i, 

- L 001 VW J "-^Ol"" ~ 1 ^02—- L J 

(2.39) a d = L d Q {Vu) d - L d l l -fcx d - L d 2 ^T 

f = L d (Viv) d - L d 21 l -f(x d - L d 2 2 -^T 

n a = Zg b (Vu;) -L8 l ^a 8 , 

(2.40) a a = LJ (Vw) a - L? 1 ^a a , 



(2.41) |a° = -L°lfi t a°. 

The pseudo scalar a and the tensors a d , and T are related to the 
distortion of the fibers, so we suppose they are of a-type, i.e., not 
changing under time inversion. On the other hand u> is of /3-type, 
i.e., changing sign under time inversion. The ONSAGER-CASIMIR 
reciprocal relations are 

(cy /io\ jd j d jd jd jd jd ja ja 

l z - 4Z J -^10 — — -^01) -^20 — ~~ -^02) -^12 — -^21) -^10 — ~ -^01- 

The inequalities 

A > 0, 7] > 0, > 0, Lf x > 0, 

(2.43) I4>0, Lg o >0, L?!>0, L° n >0 

L n L 22 > (Lf 2 ) 2 

are the consequences of the Second Law of Thermodynamics. 

According to CASIMIRs reciprocal relations, the Second Law gives 
no restrictions on the values of L d x and L 20 . The general theory does 
not imply any further restriction on the value of the coefficients but 
the analysis of some particular situations makes other simplifications 
obvious. 

Equation ()2.41|) shows that the pseudo scalar a tends to zero even if 
it were nonzero some time ago and it cannot be excited at later times, 
because there is no coupling to other variables. Therefore, after some 



Because the trace of II does not appear in the balance of energy 
and in the entropy production it is not relevant, and we conclude that 
the distinction between the symmetric and the deviatoric parts of the 
tensor II is not necessary, i.e., 

IT = U d . 

There is no coupling between the symmetric part of the velocity gra- 
dient (the affine deformation of the fluid) and the distortion of the fibers 
because d is the only second order polar tensor among the thermody- 
namic forces. If we had regarded the affine deformation of particles in 
the liquid, equation ()2.38|) should have turned to the relation known in 
visco-elasticity 

2.5. Further considerations. Up to here, the balance equations in 
general and the dissipation inequality have been used only. Now we 
consider a special flow field which is such that the distortion of 

a previously un-distorted fiber is the same as that of a fluid element, 
i.e., the equation 

(2.44) 0=Vuj-ii 

holds if ct a = ct d = T = 0. Comparing it with equation (J2.13j) and 
using equations (l2~3Hj) and (l2~4TT|) and (l2~4Tl) we get that 

Vw • n = L a 10 (Vu:) a ■ n + L d 10 (Vu>) d ■ n- 

(2.45) -L d 10 n [n ■ (Vw) rf • n] + L d 20 n [n ■ ( Vu) d ■ n] 

must hold for any n. This equation is compared with the general 
decomposition of the tensor Vu; into its symmetric part and antisym- 
metric part. As the phenomeno logical coefficients do not depend on 
the flow geometry, we conclude that 

(2.46) L a w = L d w = l and L d w = L d 20 . 

If we make the reasonable assumption that the bend relaxes with a 
single relaxation time and that the twist and the bend relax separately, 
we get that 

(2.47) 3Lf x = 5L a u and L d 12 = 0. 

These are relations between the constant (within the linear theory) 
constitutive coefficients, and therefore they are valid not only in the 
special case. 

With these results the constitutive equations (|2.39|) get the form 
n = L d (Vu) d + LS (Vu;) a - l -fcx d - l -fa a - T, 

(2.48) & = (Vu)-Lf 1 ^a 

f= (Vw) a - L d 2 ^T. 

The number of material coefficients to be determined experimentally 
has been reduced to six. 

We have shown that the stress tensor may have a skew symmetric 
part in a fiber suspension. The equations of motion are those of a 
micropolar continuum. Moreover, the theory based on the approxi- 



and the bend of the fibers relax with different relaxation times. The 
authors expect, that the improvement of the approximation (J2.13)) to 
higher orders does not change the equations because of the nature of 
the isotropic tensors. To decide if it is really so needs further investi- 
gations. 

3. Exploitation of the dissipation inequality according 

to LIU 

In this section we want to derive the restrictions on constitutive 
functions for fiber suspensions by the method of LIU JB|- From a very 
basic amendment to the second law of thermodynamics on physical ar- 
guments [63J one can show that the requirement of a positive entropy 
production density restricts possible constitutive functions and does 
not rule out certain process directions in non-equilibrium. The most 
general way to derive these restrictions on constitutive functions is the 
method by LIU ^Bj- The restrictions follow after the set of variables 
for material properties, the state space, has been chosen. They do not 
completely determine the material behavior, as it is natural, because 
there are different materials with the same set of variables, but differ- 
ent constitutive properties. The results of this method give the most 
general framework compatible with the second law of thermodynamics. 
In the example of fiber suspensions we chose the state space 

(3.1) Z={p,T,v,Vv,Q,0,fi} , 

where all constitutive quantities are assumed not to depend on v. This 
state space is related to the set of variables in the domain of the con- 
stitutive functions entropy and internal energy in the previous section. 
Apart from the equilibrium variables mass density and temperature it 
includes the velocity gradient, an orthogonal second order tensor Q 
describing the local fiber orientation (which is assumed variable here). 
Q is the mapping between the orientation vector of the (undistorted) 
fiber and a reference coordinate system. With the tensor Q it would 
also be possible to account for the orientation of biaxial elements, such 
as plates, but this is out of the scope of the present paper. It is possible 
here, that there is an arbitrary distribution of fiber orientations. Then 
the orientation distribution function can be approximated by a second 
order tensor, denoted as A. The definition of this second order tensor 
from the distribution function will be discussed in section |3J Then the 
tensor Q is the mapping between the coordinate axes and the principal 
axes of the tensor A. 

We also include the second order tensor O for the local deformation 
of the fibers. O is defined as 



(3.2) 



= \ ((VQ T ) Q) : e 



In contrast to the previous section we also include the time derivative 
Q of the fiber orientation: 

(3.3) ^ = fixQ 

in the set of relevant variables for constitutive functions. 

As the materials discussed here are micropolar media we will take 
into account the balance of internal angular momentum in addition to 
the balances of mass, momentum, and energy. The following inequality 
has to be exploited: 



+X U + V q-t:Vv-r + U: VQ^j 

(3.4) +X S ■ (p-^ - e : t + V • n T + pm J > . 



After exploiting the differentiations of the constitutive functions, de- 
fined on the state space (|3.1|) . according to the chain rule it results an 
inequality linear in the following higher derivatives: 

(3.5) f,p,v,^^,d,n,Vp,VT,VVv,VO,Vn . 
at 

These higher derivatives are not all independent, but one constraint 
between them has to be taken into account. To show this we will use 
components with respect to a cartesian coordinate system: 

Proposition: 



(3.6) '^ = { vnf + Ue.Q: dQ ^ * 



dt 2 \ dx ) dx 

or in components 



/q „, dO ik _ dVLi 1 dQ rs dv r , 



dt dxk 2 dx m dxk 



The proof of this proposition is shown in the appendix. Q, Vt> and 
VQ are state space variables, whereas Vf2 is not included in the state 



After inserting this constraint we can write down the LIU-equations, 
corresponding to the different higher derivatives: 



a,) t: P %L + ^ + P y.^o 

-> -(4-a--(- n+ i)- 

3.3) V P : » V ." y.^ + A .* „ 

ap ap ap ap 

3 . 4) vr: »_v.*-v.^ + *.*.o 

3,7) *: pA^O . 



From this set of equations the multipliers A p , A p , A s , and A" can be cal- 
culated. Between the specific spin density and the orientation change 
velocity f2 the relation 



(3.18) s = @ n 



holds. is the moment of inertia tensor, which is assumed to be 
constant. In case of suspensions of flexible fibers this is only an ap- 
proximation, because in principle this tensor changes if the fibers are 
deformed. We assume here that these deformations are small and the 
variation of © can be neglected. Then the specific spin density s de- 
pends only on the orientation change velocity f2, and all other partial 
derivatives vanish. In this case we have 



<9s 

(3.19) -- = & 



and the equations for the multipliers simplify to 
(3-20) A" = -£ = -- 



(3.23) A p = . 

By 

for the derivative -f^- we insert ■k . This relation is known in equilibrium. 

ar 

In non-equilibrium it needs some additional argumentation, which can 
be found in |64j . 

The remaining equations (I3~T!H) . (ETTTJ) . (l3~T3l . (I3~T3J) . (ETToT) . and 
(|3.1fjjl give after inserting the multipliers the restrictions on constitutive 
functions. This will be shown here only under the assumption of a 
constant moment of inertia tensor. We have: 



(3.24) 



n 




[7] — j;U) 

dVv 
d (r) — ^u) 



dO 

(3.25) : Q 

(3.26) 
(3.27) 
(3.28) 





+ 


dn 




1 dq 




dp 


T~dp 


dn 


d® 


1 dq 




df ~ 


fdf ~ 


dn 


1 


dq 






dVv TdVv dCl d\7v 



r „ 9Q , d*_idq_ d( v - i M ) an 
1 ' do Tdo~ dn do 

Introducing the free energy density f := u — Tt] and the difference 
k = <fr — the resulting restrictions on the constitutive functions can 
be written as: 

(3.30) 4t~ = 
v ' dVv 

TT _ 1 df ld<f> 1 df . dll 

(3.31) H-—P— + -— -——•(=) • — 

= fOT -^>™> • 

Equation (|3.32|) shows that the difference k = <fr — is surely 
non-zero, if the free energy density depends on the orientation change 
(M- -4- m 



variables p,T, Vv, O. In this case the very frequently made constitu- 
tive assumption (see for instance Irreversible Thermodynamics) of the 
entropy flux $ being heat flux divided by temperature is not fulfilled. 

Equation (|3.31|) is a differential equation for the couple stress. It 
reduces to an algebraic equation, if the free energy density does not 
depend on f2. In this case the couple stresses can be calculated as a 
derivative of the free energy density 



1 df 

(3.33) n = -,-L 



3.1. The entropy production. The residual inequality is built up 
by all terms in equation ([3.4)1 . which contain no higher derivatives, 
but only state space functions. The constraint equation ()3.7|) has been 
inserted, and z — ^. We will deal here only with the case of a constant 
moment of inertia tensor. For the multipliers equations (J3.20)) to (|3.23j) 
have to be inserted: 



dr] ■ p drj T d$ dQ 

a = P dQ ' Q ~ 2dO ' ' ^ ® Q "' E dQ ' ~dx 

xs fdn T dQ 

+ XU (p^--Q~^-^ v - (VQ) • Q T ■ e + (e : t) ■ s 



dQ ^ 2dO 



g.:g- t :(V., +rf .. 



1 ( du ■ p du _ / x 

dq dQ ,„ , „ \ 




change of orientational order 
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coupling between viscous flow and orientational order 





transport of orientation viscous flow 
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coupling between change of orientation and gradient of orientation 



The interpretation of the different terms is given in the equation. 

The first term in the entropy production, due to change of the ori- 
entational order is analogous to the terms involving a , a d , and T 
in the first section. Both are due to changes of the internal variable. 
The next term is new, because in this part here we have allowed for 
a dependence of the free energy density on the spatial gradient of the 
internal variable. A term due to transport of orientation does not show 
up in the entropy production according to Irreversible Thermodynam- 
ics because there the vector k is assumed to be always zero. The next 
term due to viscous flow and the term t : Vv appear in both method of 
exploiting the dissipation inequality. Terms involving the derivative of 
the free energy density with respect to the orientation change velocity 
are new, because the rotation of fibers was not taken into account in 
the Irreversible Thermodynamics treatment. On the other hand, in this 
section we were not interested in heat conduction, and the temperature 
gradient was not included in the state space. 



Fiber suspensions are an example of so called complex materials, 
meaning materials with an internal structure, which can change under 
the action of external fields resulting in complex material behavior. 
The element of internal structure in our example is the orientation and 
deformation of the fiber. There are two principally different possibili- 
ties to deal with complex materials within continuum thermodynamics: 
The first way is to introduce additional fields depending on position and 
time. These fields can be internal variables jHSEHl, order or damage 
parameters jHZ] , Cosserat-triads [UHj , directors EDI , alignment and 
conformation tensors [7_U [72] . The other way is a so called mesoscopic 



(3.34) 




4. Mesoscopic theory of fiber suspensions 



new mesoscopic fields are defined on the space M 3 x R t x M. The 
manifold M is given by the set of values the internal degree of free- 
dom can take. Therefore the choice of M depends on the complex 
material under consideration. We will see later that the manifold M 
should be such that differentiation and integration is possible on it. We 
call this way of dealing with the internal structure of complex materi- 
als a mesoscopic concept, because it includes more information than a 
macroscopic theory on x 1 () but the molecular level is not consid- 
ered like in a microscopic approach. The mesoscopic level is between 
the microscopic and the macroscopic level. The domain of the meso- 
scopic field quantities x l t x M is called mesoscopic space. The 
orientation of an undeformed fiber is described by a unit vector n, 
where turning around the fiber by tt does not change the orientation 
and therefore n — > — n is a symmetry transformation. The vector n is 
an element of the unit sphere S 2 . The deformation of the fiber is given 
by the vector tp introduced in the first part. The difference between 
that first part of the paper and the mesoscopic theory is, that now fiber 
orientation and deformation are not introduced as macroscopic quanti- 
ties, but as variables in the domain of fields. To distinguish these two 
kinds of quantities all mesoscopic quantities are denoted with a . We 
will see later that this way we can deal with different fiber orientations 
and deformations within one volume element of continuum theory. In 
our example of fibers the manifold M is the product of the set of values 
the orientation and deformation of a fiber can take: M = S 2 xR 3 with 
fi G S 2 and <p G M 3 . In order to be treated within a mesoscopic theory 
the length of the fibers must be smaller than the linear dimension of 
the continuum element, meaning that the fibers must be smaller than 
any macroscopically interesting length scale. Therefore this mesoscopic 
concept cannot be applied to fiber composite materials with long fibers. 

Beyond the use of additional variables m the mesoscopic concept 
introduces a statistical element, the so-called mesoscopic distribution 
function (MDF) generated by the different values of the mesoscopic 
variable in a volume element. Here this is a distribution of fiber orien- 
tations and deformations /(n, <p, x, t) in the volume element at position 
x and time t. This distribution function gives the probability density of 
finding a fiber of a specified orientation and deformation in this volume 
element. It is normalized 



4.1. Mesoscopic balance equations. Now fields as mass density, 
momentum density, etc. are defined on the mesoscopic space. For dis- 
tinguishing these fields from the macroscopic ones we add the word 
"mesoscopic" . For instance the mesoscopic mass density g(h, 0, x, t) 
is the mass density taking into account only fibers of a specified ori- 
entation n and deformation (p. The macroscopic mass density is the 
integral over all possible orientations and deformations: 





(4.1) 




From this equation and the interpretation of the distribution function 
as probability density it is clear that the MDF is given by the mass 
fraction: 

(4.2) f(n,<p,x,t)= qM . 

For any kind of set of mesoscopic variables balance equations can 
be derived for the mesoscopic field quantities (see for instance [23] ) 
starting out from the macroscopic global ones [221 Ell E3 Ell • A gener- 
alized Reynolds transport theorem in the mesoscopic space [73] is used 
to transform the time derivative, and Gauss theorem is applied. For 
the fiber suspensions in regular points of the continuum there result 
the local mesoscopic balance equations. We will show here only the 
balance of mass and the balance of momentum with the abbreviation 
(•) = (n, 0, x, t). Similarly mesoscopic balance equations of energy and 
of angular momentum can be derived. There is a balance of angu- 
lar momentum independently from the balance of momentum, because 
rotations of the fibers result in an internal angular momentum. 

Mesoscopic balance of mass 

d A 

w 

£>(•) is the material velocity of fibers of a specified orientation and 
length, n is the orientation change velocity, and <£?(•) is the deformation 
change velocity. 

Mesoscopic balance of momentum 

| (g(.)v(.))+V x .(v(.)g(.)v(.) - t(-)) + 

+V n • t)g(-)v(-) - f (■)) + (#(•)#(•) " *(•)) = 

(4-4) $(•)/(•)■ 

Here /(•) is the external acceleration density, £(•) the stress tensor, 
and T(-) the stress tensor on orientation space (non-convective mo- 
mentum flux in orientation space), f(-) is the momentum flux vector 
with respect to the fiber deformation variable, all quantities defined on 
the mesoscopic set of variables. 



(4.3) -£(•) + V,.(e(.)«(O)+V„-(ff(OA(O)+V # (ff(.)0(- 



0. 



As always in continuum theory the balance equations do not form 
a closed set of differential equations, but constitutive equations are 
needed, now on the mesoscopic level. We are interested here mainly in 
the constitutive equation for the stress tensor. 

4.2. Orientational order parameter and deformation variable. 

The aim is to introduce macroscopic quantities from this mesoscopic 
background, which describe the distribution of fiber orientations and 
the average distortion of fibers. These are internal variables in the 



Orientational order parameters 

(4.5) A k = / f(<p,n,x,t)n...nd 3 Lpd 2 n , 

Js 2 Jr 3 " * 

deformation order parameters: 

(4.6) $ fc =/ f f((p,n,x,t)(fi...(pd 3 ^d 2 n . 

JS 2 Jr 3 V v ' 

k 

These order parameters are tensors of successive order. They are 
macroscopic fields depending on position and time. With respect to 
fiber orientations we have the symmetry transformation n — > —n. 
Therefore all odd order orientational order parameters vanish, and the 
first non-zero order parameter, apart from the isotropic part A = 1 
is the second order one: A 2 . This can be included as an internal vari- 
able in a macroscopic exploitation of the dissipation inequality, like the 
tensor internal variable Q in section EJ 

With the mesoscopic background it is possible to derive also equa- 
tions of motion for these variables, which will be left for a future work. 

4.3. Mesoscopic and macroscopic stress tensor. The extensive 
quantity momentum has to be the integral over all mesoscopic mo- 
menta: 



(4.7) p(x,t)v(x,t) = / / p(x,t,cp ,n)v(x,t,<p ,h)d tpd n 

Js 2 Jr 3 

Integrating equation (|4.4|) over all fiber orientations and fiber defor- 
mations we must obtain the macroscopic balance of momentum (|2.2(J|) . 
This gives a relation between the mesoscopic and the macroscopic stress 
tensor. With the abbreviation (•) = (x,t,cp , n) we have 

/ / f4 (»(•)«(•)) + v, •(<>(.)»(•)*(•) -*(•)) + 

JS 2 JR3 \ at 

+V n • (ii(x,t)e(-)v(-) - f (•)) + - f (•)) d 3 0d 2 h 



(4.8) / / g(-)f(-)d^d 2 n . 

Js 2 Jrz 

Integrating the divergence terms over the whole mesoscopic space gives, 
according to GAUSS-theorem boundary terms. The orientational part 
of the mesoscopic space, the unit sphere, is a closed surface, i.e. with- 
out boundary. Concerning the deformation variable we suppose that 
deformation cannot have arbitrarily large values, and therefore there 
is no flux over the boundary at infinity. We conclude that the last two 
integrals on the left hand side vanish, and we end up with the equation 

| / / (Q(-)v(-))d 3 <pd 2 h 
at Js 2 </r 3 

+V X - f f (v(-)g(-)v(-) -£(■)) d?(pd 2 ri = 
Js 2 Jr 3 

(4.9) f [ g(-)f(-)d^d 2 h . 



Comparing the divergence terms under the space derivative with 
those in the macroscopic balance of momentum we have: 



(4.10) 



/ / («(■)£(■)£(■) -£(■)) d 3 (pd 2 h = pvv - t 
Js 2 Jm 3 



We conclude that for the constitutive quantity stress tensor the follow- 
ing relation holds 



2 h 



(4.11) t(x,t)= [ [ (£(■)- g(-)Sv(-)6v(-)) d 3 0d 

Js 2 Jm. 3 

with the abbreviation 

(4.12) 5v(-) = «(•) -v(x,t) 

This result shows that the fluxes, (here the stress tensor) are in 
general not simply the integrals of the corresponding mesoscopic quan- 
tities over the mesoscopic variables. This is true only for the extensive 
quantities. 

If all fibers have the same translational velocity, then Sv(-) = and 
the macroscopic stress tensor is the integral over all mesoscopic ones: 

(4.13) t= I [ t(-)d 3 pd 2 h . 

Js 2 JR3 

The mesoscopic stress tensor is a constitutive quantity, defined on a 
suitable set of variables. This set of variables can include mesoscopic 
quantities as well as macroscopic quantities. A reasonable assumption 
for the set of variables is: 

(4.14) Z = {p,T,n,<p,d,Vxv} . 

With this set of variables a representation theorem to linear order in 
the velocity gradient and the deformation variable (p gives the following 
expression for the mesoscopic stress tensor: 

*■ P ( 

t — - «iM + a 2 ny> + a 3 (pn + a 4 n(V x v) + a 5 CV x v)n + a 6 d 
P v 

(4.15) +a 7 nn ■ d + a%h ■ dh + agfi ■ d ■ fihn 

where the material coefficients a± to ag may all depend on the (macro- 
scopic) mass density p and temperature T. 

Averaging this over the mesoscopic variables we obtain according to 
equation ([4.13)1 : 

t= / — (oiihh + a 2 n<£> + a 3 ipn + a 4 n(V x v) + a 5 (V x v)h 
Js 2 Jm 3 P 

+a§d + ajfin ■ d + a&ri ■ dh + agfi ■ d ■ hhfi) d z (pd 2 n . 

= aiA 2 + a 2 (ncp) + a 3 (cpn) + a e dd 3 (pd 2 n 



The average of «4n(V x v) vanishes, because J s2 fnd 2 n = due to the 
symmetry n <-> — n, analogously for the term with a$. The averages 
(flip) and (0n) are non-zero, because they are even functions of h: 

(4.17) 0(-fi) = -0(h) 
(see the first section) and therefore 

(4.18) — h(f(— h) = h<p(h) . 

The stress tensor (|4.1fi|) clearly can have an antisymmetric part t a : 
(4.19) 

1 1 / = o \ 

t a = - (a 2 - a 3 ) ((hep) - (ipn)) + -(a 7 — a s ) (A - d- d - Ay . 

5. Conclusions 

We have investigated the constitutive properties of suspensions of 
long fibers, which can be deformed. In the second and in the third 
part we have taken into consideration in addition the possibility of dif- 
ferent fiber orientations. This complex material has been treated with 
three different methods. Two of these methods are different ways of 
exploiting the second law of thermodynamics in macroscopic contin- 
uum thermodynamics. The first method was classical Thermodynam- 
ics of Irreversible Processes (TIP). The second method of exploitation, 
the method according to LIU, is more related to Rational Thermody- 
namics. In both methods the aim is a derivation of information on 
constitutive functions, and the first method shows that in fiber suspen- 
sions there can exist an antisymmetric part of the stress tensor due to 
internal degrees of freedom. The method of TIP gives additional relax- 
ation equations for the internal variables introduced in the beginning 
in order to account for the internal structure. Such relaxation equa- 
tions cannot be derived by the method of LIU. The assumptions made 
in the beginning by this method are less restrictive, and the constitu- 
tive equations consistent with the results of the method by LIU are 
more general than the equations derived from TIP. It could be shown 
that the assumption made in TIP concerning the entropy flux being 
heat flux divided by temperature does not contradict the results of the 
LIU-procedure, if the free energy density does not depend on the ori- 
entation change velocity of the fibers. This assumption was made by 
the choice of variables in the first section. Therefore this exploitation 
is in agreement with the results of the LIU-procedure. However, if 
the orientation change velocity is included in the set of variables, the 
entropy flux might be not simply heat flux divided by temperature, 
as a result of the exploitation of the dissipation inequality. This can 
be interpreted as non-convective entropy transport due to change of 
fiber orientations. In addition couple stress could be calculated as a 
derivative of the free energy. 

Finally the idea of a so called mesoscopic theory was sketched. There 
we introduced a finer description considering single fibers with differ- 
ent orientations and different states of deformation. The macroscopic 



over the different fiber orientations and different states of deformation. 
This method gives, apart from the usual balance equations, a definition 
of internal variables in terms of averages. Equations of motion for these 
internal variables will be discussed in a future paper. Another result 
shown here is the fact that the stress tensor calculated as an average 
of mesoscopic quantities has an antisymmetric part. 

In the mesoscopic theory we did not deal further with the possibil- 
ity that the orientation distribution of fibers and also the distribution 
of fiber deformations can change, for instance under the action of a 
flow field. It is expected that both effects have influence on material 
properties, a problem, that will be dealt with in a future work. 
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Appendix: Proof of the proposition on the time 
derivative of the state space variable O 



Proposition: 



d ° /„^t i ( „ d Q T \ dv 
(5.1) — = (vny+-(e-Q 



dt 2 V dx J dx 

or in components 

/_ x dO ik dVti 1 dQ rs dv r , 

(5.2) — - = — + -e irl Qi s - 



dt dxk 2 dx m dxk 



Proof: 

In components we have 
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